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> ■ ABSTRACT 

£*""». ' We have compiled a sample of 2728 nearby (z < 0.08) elliptical galaxies with pho- 

tometry in the g, r, i, z bands from the Sloan Digital Sky Survey (SDSS) and J, H, K 
photometry from the Two Micron All Sky Survey (2MASS). Stellar masses, stellar 
velocity dispersions and structural parameters such as sizes and surface mass densi- 
ties are also available for these objects. In order to correct the aperture mismatch 
between SDSS and 2MASS, we correct the SDSS magnitudes to the isophotal circu- 
lar radius where the 2MASS magnitudes are measured. We compare the correlations 
• between optical, optical-infrared and infrared colours and galaxy luminosity, stellar 

Q-f mass, velocity dispersion and surface mass density. We find that all galaxy colours 

correlate more strongly with stellar mass and velocity dispersion than with any other 
£h ' structural parameter. The dispersion about these two relations is also smaller. We also 

^ \ study the correlations between a variety of stellar absorption line indices and the same 

set of galaxy parameters and we reach very similar conclusions. Finally, we analyze 
correlations between absorption line indices and colour. Our results suggest that the 
optical colours of elliptical galaxies are sensitive to a combination of age, metallicity 
■ and a-enhancement, while the optical-infrared colours are sensitive to metallicity and 

5— i ' to a-enhancement, but are somewhat less sensitive to age. 
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1 INTRODUCTION 

It is well known that the observed properties of el- 
liptical galaxies obey a number of tight relations: the 
most well-studied of these are the colour-magnitude rela- 
tion (CM R), the fundamental plane and the Mg-cr rela- 
tion <|Faber fc Jackson] 1976|:|Visvanathan fe Sandage|l977[ 
[Dressier et al.ll 1987tlDiorgovski fc Davislll987j:l Bower et alJ 
Il992l: lEllis et alJll997t iBernardi et alJ l2003alibl Jdft . These 
correlations link the properties of the stellar populations of 
early-type galaxies with their masses and their structural 
properties and are believed to encode important informa- 
tion about how early-type galaxies formed. 

It has long been known that early-type galaxies in 
nearby clusters exhibit a tight colour-magnitude relation, i.e. 
more luminous ellipticals have redder colours than less lumi- 
nous ellipticals l|Visvanathan fc Sandagelfl977t iBower et alJ 
1992). This work has been extended to moderate red- 
shift and it is found that there is no significant change 



in the slope with redshift lAra^c4^alam^nca^^L|_^9d 
[ Stanford etallll995l. ll99Sl : lEllis et al.lll997l : iKodama etaH 
\199& iHolden et alJ l20oll . The zero-point of the relation 
shows modest colour evolution, consistent with the passive 
ageing of a stellar population formed before z ~ 2. Further- 
more, the slope of the CMR of ellipticals has been found 
to show little dependence on environment ([Terlevich et alJ 
1200 It IBernardi et aljl2003dl: iHogg et"afll2004 . 

The origin of the CMR is still a matter of contro- 
versy. The conventional interpretatio n of the CMR is that 
it is primarily a metall i city effect |Faber | [l977|: |DressleJ 
[ 19841: lArimoto fc Yoshill Il987t IKodama fc Arimotol I199tF 
lFerrera^^t^Ln^99F ~R is often assumed that elliptical 
galaxies form monolithically in a single giant burst of 
star formation at high redshift. Since the binding en- 
ergy per unit mass of gas is higher in more massive 
galaxies, they are able to retain their gas for longer and 
so reach higher metallicities than less massive galaxies 
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<LarsorJll974l:lArimoto fc Yoshiilll987l: iBressan et alJ|l99d 
iKodama fc Arimotol |1997T) . IKauffmann fc Chariot! (Il998l) 
and iDe Tucia^^Ll(|20o3 ) have argued that the CMR can 
also be explained in hierarchical models of galaxy formation 
when strong feedback is included. In these models, the CMR 
for elliptical galaxies arises because large ellipticals form by 
mergers of massive, metal-rich progenitor disk systems. In- 



deed, some evidence of that buil d-up in the early-type g alax 
population has been observed llChen et aljl2p03l iBell et~a 
I 2004 ICross et all 12004 iFerreras et alJl2005ft . However, i 
is also possible to explain at least part of the CMR as 



an increase in mean stella r age with lu minosity l|W or^hevl 
119961: IFerreras et al]ll999l) . Furthermore, lj(3rgenseiT^^99Tl 
and iTraeer et al.l (|2p00) find strong evidence for a signifi- 
cant intermediate-age population in some elliptical galaxies. 

Additional information about stellar populations can 
be gained by studying the strength of stellar absorption 
features in galaxy spectra. The standard way of quantify- 
ing absorption line strengths in galaxies is via the Lick-ID S 
system, which includes 25 spectral absorption features de- 
fined by a 'feature' bandpass and two 'pseudo-continuum' 
bandpasses. The indices are calibrated using the spectra of 
460 Galactic stars over the wavelength range from 4000A to 
6400 A. The sensitivities of different Lick indices to stellar 
age a nd metallicity have been discussed in a number of pa- 
ers (Burstein ct al 1ll984fc iGoreas et alJll99i iTraeer et alJ 



ll99Sl:lHenrv fc Worthevlll999l) . The H-Balmer lines are pri- 
marily age-sensitive. Likewise, the 4000A break (i.e. the ra- 
tio of the flux blueward and redward of 4000A) gets stronger 
with age, but for old stellar populations, it shows also a 
second-order sensitivity to metallicity. Indices which are pri- 
marily sensitive to metallicity include Fe and Mg lines be- 
tween 4500A and 5700A and molecular features such as 
the TiO band. In massive elliptical galaxies, the abundance 
of a-elements with respect to Fe-pe ak elements can diffe r 
from the scaled-solar abundance ratio dWorthev et al.ll992>) . 
Indices which trace the abundance of either a-elements or 
Fe-peak element s can therefore be se nsitive to the degree of 
a-enhancement jThomas et al.ll2003l) . 

In addition to correlations with galaxy luminosity, past 
studies of elliptical galaxies have also focused on how colours 
and absorption line indices depend on the stellar velocity 
dispersion of the galaxy. To first order, both the luminos- 
ity and the velocity dispersion of an elliptical are measures 
of its mass; the luminosity measures the mass contained in 
stars and the velocity dispersion is a dynamical probe of 
the total mass contained within some radius. The star for- 
mation history of a galaxy, however, m a y not depend pri- 
marily on its mass. Indeed, iKennicuttl l)l998l) has shown 
that star formation rates in spiral galaxies correlate best 
with the local surface density of gas and a law of the form 



1Mqpc~ 



'MQyr 1 kpc 2 pro- 



Esfr = (2.5 ±0.7) X (_ 

vides a good fit to the data. IBell fc de Jond <|2000T) stud- 
ied trends in the optical and infrared colours of a sample 
of nearby spiral galaxies and found that the colours most 
sensitive to star formation histor y correlate best w i th stel- 
lar surface density. Furthermore. IKauffmann et al] (|2003b) 
studied the relations between stellar mass, star formation 
history, size and internal structure of a complete sample of 
10 5 galaxies and found that the fi, — M* relation for the late- 
type galaxies can be well described by a single power law. 



They also found that the star formation history of low-mass 
galaxies are correlated more strongly with surface density 
than stellar mass. These results suggest that local rather 
than global factors regulate the rate at which spiral galaxies 
form their stars at the present day. Motivated by the sur- 
face density dependence of spiral star formation histories, it 
seems worthwhile to search for surface density dependence 
in the properties of early-ty pe galaxies. 

Finally, in a recent work lEisenstein et al] i2003l) studied 
how the average optical spectra of massive galaxies in the 
SDSS vary as a function of both luminosity and environ- 
ment. They found that as redshift, luminosity and environ- 
ment change, the variation of pairs of Lick Indices appears to 
follow a simple one-dimensional locus, suggesting that vari- 
ations in these three parameters do not lead to independent 
changes in the spectra of elliptical galaxies. 

In this paper, we present the correlations between the 
colours of elliptical galaxies and a range of different struc- 
tural parameters. We make use of data from the Sloan Dig- 
ital Sky Survey (SDSS) and the Two Micron All Sky Sur- 
vey (2MASS) to construct a sample of early-type galaxies 
with photometry in five optical (u, g, r, i, z) and three near- 
infrared bands J (1.25/im), H (1.65/^m), K 3 (2.17/im) bands. 
Stellar velocity dispersions and size measurements for these 
systems are available from the SDSS datab ase and estimates 
of stella r mass are taken from the work of IKauffmann et al] 
J2003al) . The outline of the paper is as follows. Section 2 
describes how our sample is selected, the method adopted 
for doing the aperture correction and calculating (K+E)- 
correction of magnitudes. Section 3.1 presents the correla- 
tions between optical, optical-infrared, infrared colours and 
structural parameters. In section 3.2, we study the relations 
between colours and a variety of Lick indices. Finally, in sec- 
tion 4 we summarize our findings and give our conclusions. 



2 THE SAMPLE 
2.1 Sample selection 

The SDSS is an imaging and spectroscopic survey of the 
high Galactic latitude sky visible from the Northern hemi- 
sphere, which will obtain u, g, r, i and z photometry of al- 
most a quarter of the sky and spectra of at least 700,000 
objects jYork et al] 1200(11 . A description of the soft ware 
and data produc t s are g iven i n Istoughton et al] J2002T) and 
lAbazaiian et al] feOOl l2004tl . The 2MASS is a ground- 
based, near-infrared imaging survey of the whole sky and 
its extended source catalog (XSC ) contains almost 1.6 mil- 
lion galaxies Jjarrett et aljboOOallrl . We positionally match 
sources observed by the 2MASS to galaxies in the 'main' 
spectroscopic sample of the SDSS Data Release One (DR1; 
Abazajian et al. 2003). Galaxies in this sample have 14.5 < 
r < 17.77. Practically all 2MASS sources in the DR1 area 
(~ 97% for extended sources from the XSC) are matched to 
an SDSS counterpart within 2 arcsec. We restrict our anal- 
ysis to a narrow range of redshift 0.02 < z < 0.08 in order 
to minimize uncertainties introduced by K-corrections. The 
entire matched sample contains 8166 galaxies in this range 
of redshift. 

In this paper, we focus on early-type galaxies. To se- 
lect these systems, we require that they have 4000 A break 
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strength greater t han 1.6 (we adopt the narrow definition of 
the break given in iBaloeh et al .119991) and concentration pa- 
rameter C = i?90,r/-R50,r > 2.6 (where Rgo,r and Rso,r are 
the radii enclos ing 90% and 50% of the to tal Petrosian light 
of the galaxy). iKauffmann et aD i2003tJ) have shown that 
this cut separates massive, dense, early-type galaxies from 
low mass, low density late-type systems. In order to exclude 
systems with ongoing star formation or residual AGN activ- 
ity, we also eliminate all galaxies where the equivalent width 
of the Ha emission line is greater than 2 A . Our final sample 
contains 2728 early-type galaxies. 



2.2 Photometric quantities and aperture 
corrections 

In order to compare optical and optical-infrared colours, it 
is important that both SDSS and 2MASS magnitudes are 
measured within the same aperture. There are several sets of 
magnitudes provided in the 2MASS XSC: elliptical isopho- 
tal magnitudes, circular isophotal magnitudes, total mag- 
nitudes using a Kron elliptical aperture, total magnitudes 
from extrapolating the fit to the radial profile and magni- 
tudes within fixed circular apertures (5, 7, 10, 15 arcsec ). In 
this paper, we use the isophotal fiducial magnitudes, which 
are measured within the circular aperture corresponding to 
a surface brightness of 20.0 mag arcsec -2 in the Jf s -band 
(the aperture is denoted as Rk20fc)- At this isophotal mag- 
nitude, the background noise is still relatively low, but the 
aperture is large enough to enclose most of the light from 
the galaxy (see the 4th panel of Figure 1). This set of magni- 
tudes are measured with the same size aperture in all three 
(J, H, K a ) 2MASS bands and they thus provide a consistent 
set of infrared colours. 

SDSS provides the azimuthally averaged surface bright- 
ness in a series of circular annuli. In order to derive optical- 
infrared colours, we have calculated g, r, i, z magnitudes that 
are matched to the 2MASS measurements by interpolating 
the cumulative radial surface brightness profile in each SDSS 
band at the corresponding isophotal radius Rk20fc- From 
now on, we refer to these corrected magnitudes as the "SDSS 
magnitudes" and they are used to estimate the optical and 
optical-infrared colours throughout the paper. 

In Fig. Q we plot the distribution of the r-band Pet- 
rosian radius (from SDSS), the isophotal radius Rk20fc (from 
2MASS), the quantity 2Rp — Rk20fc, and the r-band flux 
within Rk20fc divided by the Petrosian flux. A small fraction 
of galaxies (about 1.2 percent) have isophotal radius larger 
than 2Rp. We simply adopt the Petrosian magnitudes for 
these galaxies and make no aperture correction. Fig.Qshows 
that for most galaxies, a large fraction of Petrosian flux is 
contained within Rk20fc - Because colour gradients are small 
for early-type galaxies l)Wu et al]l200Rh . it is safe to treat 
the colours deduced from the corrected SDSS magnitudes 
as a global galaxy colours. 

Stellar mass estimates for t his sam ple are available from 
the work of K auffmann et al ] J2003al) . To estimate stellar 
masses, two stellar absorption line indices, the 4000 A break 
strength D„(4000) and the Balmer absorption line index 
US a, are used to constrain the mean stellar age of a galaxy 
and the fraction of its stellar mass formed in bursts over 
the past few Gyr. A comparison with broad band magni- 
tudes then yields estimates of dust attenuation and of stellar 
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Figure 1. The distributions of the Petrosian radius in r-band 
R P (from SDSS), the isophotal radius Rk20fc (from 2MASS), 
the quantity 2Rp — Rk20fc an d the r-band flux within Rk20fc 
divided by the Petrosian flux (i.e., within 2Rp). 



mass. The optical photometric quantities, such as Petrosian 
magnitudes in 5 bands, size Rso,z, concentration parameter 
C = Rgo,r/R50,r, ar e obtained from the F irst Data Release 
(DR1) of the SDSS llAbazaiian et alJl2005D . Velocity disper- 
sion measurements for each galaxy are also available from 
the database and are described in Schlegel et al (in prepa- 
ration). 

Throughout the paper, we assume a cosmology with 
M = 0.3, Q A = 0.7 and H =70 km s" 1 Mpc" 1 . 



2.3 (K+E)-correction 

When converting from the apparent magnitude to the rest- 
frame absolute magnitude of an object, the fact that the 
SDSS and 2MASS filters measure the light from a fixed spec- 
tral range in the local frame needs to be taken into account. 
Correcting for this effect is known as the K-correction. The 
E-correction accounts for the change in the galaxy's lumi- 
nosity between the time that its light was emitted and the 
present day. In this paper, the sum of K-correction and E- 
correction is denoted as the (K+E)-correction. 

In order to evaluate these corrections, we generate a 
grid of model spectral energy dist ributions (SEDs) using the 
population synthesis models of jBruzual fc Chariot] Eool. 
BC03). We ado pt the universal initial mass function from 
IChabrie'rl i2003fl and lower and upper mass cut-offs are taken 
as rtiL=O.lM0 and mfy = 100Mo, respectively. We parame- 
terize the models according to a metallicity Z and the star 
formation timescale r. Galaxies form stars according to the 
equation 



(1) 



for a total duration of 12 Gyr. The effect of extinction on the 
SEDs is not included, because the d ust attenuation measure- 
ments of [Kauffmann et alJ l|2003ali show that there is very 
little dust in these systems. For each model, we store ob- 
server frame colours at a series of closely spaced redshifts 
from z=0.08 to z=0.02. For each galaxy in our sample, we 
find the model that most closely reproduces its observed 
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Table 1. The results of linear fits. The first set of quantities are the best-fit slopes. The second set of quantities disp give the dispersions 
about the relation, which we define as ^' n H ' ^ , where y° is the observed colour and j/*™ is the colour predicted by the fit and n is 

the number of galaxies. The third and fourth set of quantities AI measure the changes in colour over the interval in magnitude, mass, 
a or density that contains 90% of the galaxies (we exclude the lower and upper 5th percentiles of the distribution). AI are normalized 
by dividing either by the dispersion (set 3) or by the total range in colour enclosing 90% of the galaxies (set 4). 
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colours (including g — r,r — i,r — z,r — J,r — K and J — K) 
and we use the model to transform to colours and magni- 
tudes denned at z=0.04. Because of the very limited range 
of redshifts of the galaxies in our sample, these colour cor- 
rections are small: 0.04 mag for g — r, 0.01 mag for r — i and 
r — z, 0.03 mag for r — J, 0.05 mag for r ~ K and 0.02 mag 
for J — K on average. 



3 RESULTS 

3.1 Relations between colours and structural 
parameters 

In this section, we present the correlations of the optical, 
optical-infrared and infrared colours of the ellipticals in our 
sample with the following parameters: 

(i) absolute r-band magnitude, M r (Note that we use Pet- 
rosian magnitudes to calculate M r ); 

(ii) stellar mass M„; 

(iii) measured stellar velocity dispersion o meas \ 

(iv) "inferred" stellar velocity dispersion: <j; n / = 
C(GM*/R 50 , Z ) 1/2 , where C was chosen so that the aver- 
age value of <7i n f for all the galaxies in the sample would be 



equal to the average value of a for the same set of galaxies 
(we find C = 0.5); 

(v) stellar surface mass density /j,* = 0.5M*/(7r_R§ ,z); 

fvi) stellar mass "density" p„ = 0.5M,/(4/37r_R| 0i2 ). 

where R$o,z is the radius enclosing 50 percent of Petrosian 
flux in z-band (and hence approximately 50 percent of the 
stellar mass). 1 In this paper, the parameters M», a m eas, 
ffinf, M*> P* are m units of Mq, km/s, km/s, Mo/kpc 2 and 
M0/kpc 3 , respectively. The results are presented in Figures 
2-5, where we show the correlations between g — r, r — i, r — J 
and J — K colours as a function of these parameters. We 
find that a simple least-squares fit to the data is sensitive to 
whether or not we exclude galaxies that lie far away from the 
main relations. We thus switched to robust techniques and 
we fit a straight line to the data by minimizing the absolute 
deviations. This is shown as a solid line on the plots. 

Table 1 lists the parameters of the fits. The first set of 
quantities are the best-fit slopes. The second set of quanti- 



1 We also tried different measurements of z-band half-light radius 
including the se eing deconvolved Ser sic and Petrosian half-light 
radii derived by iBlanton et alj 1200311 . We obtained almost iden- 
tical results. 
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Figure 2. The correlations between g — r and r-band magnitude M r , stellar mass M, (in units of Mq), measured stellar velocity 
dispersion a m eas (in units of km/s), "inferred" stellar velocity dispersion a in f = 0.5(GAf*/i?5o jZ ) 1//2 (in units of km/s), stellar surface 
mass density fi, = 0.5M*/(7rii| z ) (in unit of MQ/kpc 2 ) and density p* = 0.5M*/(4/37ri?g ) (in units of Mo/kpc 3 ). The solid lines 
show linear fits using robust techniques. 



ties disp give the dispersions about the relation, which we 
define as Y17=i Vi ; wnere J/° is the observed colour and 
j/™ is the colour predicted by the fit and n is the number 
of galaxies. The third and fourth set of quantities A7 mea- 
sure the changes in colour over the interval in magnitude, 
mass, a or density that contains 90% of the galaxies (we 
exclude the lower and upper 5th percentiles of the distri- 
bution). In order to compare these colour changes between 
different photometric bands, we normalize AI by dividing 
either by the dispersion (AI /disp) or by the total range in 
colour enclosing 90% of the galaxies (AI /range). 



We find that the optical (g — r, r — i, r — z), optical- 
infrared (r — J, r — K) and infrared (J — K) colours all 
correlate with r-band absolute magnitude, i.e. more lumi- 
nous ellipticals tend to have redder colours than less lumi- 
nous ellipticals. Table 1 shows that colours correlate more 
strongly with stellar mass and velocity dispersion than with 
r-band magnitude and that the dispersions about the rela- 
ti ons are also sma l ler. Our resu lts are consistent with those 
of lBernardi et al ] J2003c]l 2005), who demonstrate that the 
colour-magnitude relation is a consequence of the fact that 
both the luminosity and the colours of early-type galaxies 
are correlated with their stellar velocity dispersions. These 
results suggest that the colours of elliptical galaxies are pri- 
marily determined by their mass. 



3.2 Stellar absorption features 

In this section we study the relations between stellar ab- 
sorption features, colours and structural parameters of the 
galaxies in our sample. 

We have measured Lick indices from the SPSS spec- 
tra, u sing the bandpass definition of IWorthev fc Ottavianj 
(1997). It should be noted that the strengths of some of 
the absorption indices are strongly affected by broadening 
due to the velocity dispersion a of the stars in the galaxy. 
Indices measured with the narrowest 'pseudo-continuum' 
bandpasses show the strongest dependence on a. To prop- 
erly compare index strengths of galaxies with different ve- 
locity dispersion and study their intrinsic dependence on 
various galaxy properties, we have to correct for this effect. 
We normalize the index strengths of all galaxies to a com- 
mon velocity dispersion of 200 km s~ x . 2 This is achieved 
by fitting a relation between index strength and velocity 
dispersion for BC03 simple stellar populations of different 
metallicities. These relations are then interpolated to the 
metallicity of the gal axy, as estimated using the methods of 
iGallazzi et all (12005ft . 

In Table 2 we summarize the results for a variety of 
different Lick indices and for g — r and r — J colours. As dis- 



2 The average velocity dispersion of the galaxies in the sample is 
160 km s _1 . 
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Figure 3. Same as figure 2, but for r — i colour. 
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Figure 4. Same as figure 2, but for r — J colour. 
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Figure 5. Same as figure 2, but for J — K colour. 




Figure 6. Spectral absorption features (corrected for velocity dispersion) of the galaxies in our sample plotted as a function of g — r 
and r — J colours. Here we show only a subset of the 29 spectral indices analyzed in the paper. The indices are ordered according to the 
significance of the correlation (given by the correlation coefficient in each panel) with g — r. 
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Table 2. Relations between index strengths (corrected for velocity dispersion) and g — r and r — J colours, fitted with our robust 
technique on the subsample of galaxies that lie within 2.5 times the dispersion in the g — r, M r and r — J, M r colour-magnitude relations. 
Column 2 gives the average S/N ratio of each index, calculated as the ratio between the index strength and the observational error. 
For each colour the first and second columns give the slope and rms scatter of the relation. The third column indicates the correlation 
coefficient according to a Spearman rank-order test. The fourth column gives the change in index strength (along the fitted relation) 
over the interval in colour that contains 90 percent of the points (A I), normalized to the dispersion about the fitted relation (disp). 
Columns f I and f 2 provide some indications about the indices. Column If lists the elements that are expected to contribute most to the 
index strength (in the case of iron, a '+' indicates a weak influence and a 'H — h' a strong influence). Column 12 indicates the sensitivity of 
the index to a/Fe variations, as expected from Henry & Worthey (1999), Thomas et al. (2003) and Thomas et al (2004). Up and down 
arrows, respectively, correspond to increasing and decreasing index strengths with increasing c?/Fe ratio, slanted arrows to intermediate 
trends, and a hyphen to no sensitivity. 
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cussed bv lThomas et all <l2003t) . the a-element to Fe abun- 
dance ratio can be empirically quantified by the ratio be- 
tween Mg and Fe indices. Here we use both the Mg2/(Fe) 
and the Mgb/(Fe) index ratio. 3 It is also possible to define 
other combinations of indices that are almost independent 
of q-enha ncement. Those inves tigated here are [MgFe] ' , as 
defined bv lThomas et alJ (120031) . and [MgiFe] and [Mg 2 Fe], 
as given by BC03. For each index, we list in Table 2 the 
slope of the relation (columns 3 and 7) , the dispersion about 
the relation (columns 4 and 8), the significance of the rela- 
tion according to a Spearman rank-order test (columns 5 
and 9) and the variation in index strength over the colour 
interval that contains 90 percent of the points (Al), normal- 
ized to the dispersion in the fitted relation (columns 6 and 



3 The (Fe) is the average of Fe5270 and Fe5335 index strengths. 



10). There is a small fraction of outliers in g — r and r — J 
colours (probably galaxies with large systematic uncertain- 
ties in their photometry) which significantly affect the fits. 
To account for this, we decide to exclude those galaxies that 
lie above 2.5 times the dispersion in the colour-magnitude 
relations shown in Figs. |5|and2] In this way, about 12 per- 
cent of the galaxies are excluded from the fit. We than adopt 
the same robust technique to fit relations between index 
strengths and colours. We consider as statistically not sig- 
nificant those relations with a correlation coefficient lower 
than 0.3. In Table 2 we also provide information about the 
indices, including the mean signal-to-noise ratio of the index 
(column 2), the elements that are expected to contribute 
most to the index strength (column 11) and the sensitivity 
of the inde x to a/Fe ratio (column 12 ). This information was 
tak en fromjHenrv fc Worthed (Il999l) : iThomas et alJ (l2003h 
and lThomas et al] J2004D . In Fig.lol we illustrate a selection 
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Table 3. Relation between index strengths (corrected for velocity dispersion) and structural parameters: stellar mass (logM*), velocity 
dispersion (log cr meas ) and surface stellar mass density (log/i,), fitted by applying our robust technique to the subsample of galaxies 
that lie within 2.5 times the dispersion in the g — r,M r and r — J,M r colour-magnitude relations. Columns 2,6,10 give the slope of 
the relation; columns 3,7,11 give the rms scatter about the fitted relation; columns 4,8,12 give the correlation coefficient according to 
a Spearman rank-order test. Columns 5,9,13 give the change in index strength along the fitted relation over the interval in structural 
parameter that contains 90 percent of the points, normalized to the dispersion about the fitted relation. 
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of (statistically significant) relations between index strength 
and g — r, r — J colours. Indices are ordered from top to bot- 
tom with decreasing correlation coefficient (indicated in each 
panel) with respect to g — r. 

In Table 2, the quantities of most interest are the 
relative strengths of the correlations between index strength 
and colour for different indices (as given by the correlation 
coefficient and, e.g., A//disp). For g — r, the strongest 
correlations are found for the Mg- and CN-features, 
the 4000A break D n (4000), the Balmer absorption lines 
(H7A, HJa and H/3) and the a-enhancement independent 
indices [Mg2Fe], [MgFe]' and [MgiFe]. When we compar- 
ing these correlations with those involving r — J we see that: 

(i) D n (4000) and the Balmer absorption lines correlate 
more strongly with the optical colour than with the optical- 
infrared colour. 

(ii) The a/Fe independent Mg-Fe composite indices cor- 
relate better with optical-infrared than with optical colour. 

(iii) Mg features and NaD correlate equally well with both 
colours. 



(iv) Mg2/(Fe) shows a weak correlation with both 
colours. 

In summary, both optical and optical-infrared colours are 
sensitive to metallicity and to element abundance ratios. 
However, the age of the stellar populations (as given by 
D n (4000) and Balmer absorption lines) influence the optical 
colours more strongly than the optical-infrared colours. One 
caveat is that the colours are measured for the galaxy as a 
whole, whereas the Lick indices are measured within the 3 
arcsecond fibre aperture. Gradients in index strength may 
thus affect the analysis. We have checked whether our results 
change if we use SDSS "fibre" colours (calculated within a 
2.5 arcsec aperture) rather than the Petrosian colours. There 
are changes in the strengths of some of the correlations, but 
our main qualitative conclusions remain unchanged. 

We now turn to the correlations between index 
strengths and structural parameters, such as stellar mass 
(logM*), measured velocity dispersion (log o- mea s) and sur- 
face mass density (log fi* ) . For consistency with the previous 
analysis, we apply here our robust fitting technique only on 
galaxies within 2.5 x disp of the g — r, M r and r — J, M r re- 
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Figure 7. Spectral absorption features (corrected for velocity dispersion) of the galaxies in our sample plotted as a function of the 
measured velocity dispersion, the stellar mass and the stellar surface mass density. The same subset of spectral features of Fig.HJis shown 
here. The indices are ordered according to the significance of the correlation with velocity dispersion. In the middle and right panels, the 
dot-dashed lines show the linear relations for galaxies with similar velocity dispersion (cyan: log<r mea3 < 2, blue: 2 log<r mea3 < 2.1, 
green: 2.1 < log a meas < 2.2, red: 2.2 < log <r meas < 2.3, magenta: logowas ^ 2.3). 



lations. The results of the fits for all the indices are summa- 
rized in Table 3, where we give slope, dispersion, correlation 
coefficient and A//disp as before. Fig.Qshows a selection of 
absorption indices against velocity dispersion, stellar mass 
and surface mass density, ordered with decreasing correla- 
tion coefficient with respect to velocity dispersion. 

The main result of this analysis is that all indices corre- 
late much more strongly with velocity dispersion than with 
any other structural parameter. This is another strong indi- 
cation that the properties of elliptical galaxies depend most 
strongly on the velocity dispersion of the system. We note 
that both the velocity dispersion and the spectral indices 
are measured within the fibre, while stellar mass and sur- 
face mass density pertain to the galaxy as a whole. One 
might thus worry that this might artificially strengthen the 
correlation between stellar absorption indices and velocity 
dispersion. We have checked that similar results are obtained 
when we correlate the absorption indices with the inferred 
velocity dispersion obtained by diving the stellar mass by 
the Petrosian half-light radius. We conclude that our main 
result is robust against these effects. 

Finally, we have checked whether there are residual cor- 
relations between age-dependent and metallicity-dependent 
spectral indices and structural parameters such as surface 
mass density. The dot-dashed lines in the middle and right 
panels of Fig. |7| show the linear relations with stellar mass 
and surface stellar mass density, for galaxies in five bins of 
velocity dispersion. This clearly shows that, after removing 
the dependence on velocity dispersion, the absorption line 



strengths do not show any significant residual correlation 
with stellar mass and, even more so, with sur face stellar 
mass density. Kennicutt's law of star formation ( Ke nnicuttl 
1998) relates the star formation rate in galactic disks to the 
surface density of cold gas in the galaxy. The fact that we see 
very little correlation between spectral indices and surface 
density for ellipticals suggests that global processes, related 
to the mass of the system, have been more instrumental in 
setting the properties of galactic spheroids. 



4 DISCUSSION 

The main aim of this paper is to investigate what can 
be learned from the optical, optical-infrared and infrared 
colours of elliptical galaxies. We have put together a sam- 
ple of ellipticals from a matched sample of galaxies with 
photometry from both SDSS and 2MASS and spectroscopy 
from SDSS. We study the correlations between colours and 
a wide variety of different structural parameters, as well as 
the correlations between stellar absorption features, colours 
and structural parameters. We find that: 

(i) Luminous and massive elliptical galaxies have redder 
optical, optical-infrared and infrared colours than less lumi- 
nous and lower mass ellipticals. 

(ii) The optical, optical-infrared and infrared colours of 
elliptical galaxies correlate more strongly with stellar mass 
and velocity dispersion than with other structural parame- 
ters. The dispersion about these two relations is also smaller. 
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(iii) Both optical and optical-infrared colours are sensitive 
to metallicity and to element abundance ratios. However, the 
age of the stellar populations influences the optical colours 
somewhat more strongly than the optical-infrared colours. 

Moreover, our results show that all Lick indices corre- 
late more strongly with velocity dispersion, either as mea- 
sured within the fibre or as quantified by the ratio between 
stellar mass and Petrosian half-light radius, than with any 
other stru c tural parameter studied. As demonstrated by 
ICappellarH j2005l) . the velocity dispersion scaled by the size 
of the system (in this case, the effective radius) provides a 
remarkably good estimate of the total (stellar plus dark mat- 
ter) mass of the galaxy. We have also studied correlations 
between colours, spectral indices and the mass inferred from 
the measured velocity dispersion Mtot = (2(j meas ) 2 i?5o, z /G. 
In Fig. [H]we compare these relations with those against stel- 
lar mass for g — r colour, a representative age-sensitive in- 
dex (Hja) and a representative metallicity-sensitive index 
(Mg2/(Fe)). In the left-hand panels the dot-dashed lines 
show the relations with total mass at fixed stellar mass, 
while in the right-hand panels they show the relations with 
stellar mass at fixed total mass. The middle and bottom 
panels show that the relation with total mass for galaxies 
with similar stellar mass has almost the same slope as the 
relation obtained for the sample as a whole. On the contrary, 
at fixed total mass, the relation with stellar mass becomes 
much shallower. 

However, the dependence on these two physical param- 
eters is less distinguishable in the case of g — r colour. It is 
difficult to tell from Fig. |H| which is the most fundamental 
parameter between stellar and total mass. Since both the ab- 
sorption features and the velocity dispersion (which enters 
in the estimate of total mass) are measured within the fiber, 
while stellar mass is a global quantity, one might worry that 
the relation with total mass is strengthened by aperture ef- 
fects. We checked the correlation between colours, spectral 
features and stellar mass in the fiber M t jiter (obtained by 
scaling stellar mass M* with the ratio between fiber and 
Petrosian z-band luminosity) and show part of the results 
in Table 01 It is shown that both the colours and spectral 
indices correlate more strongly with the mass inferred from 
the velocity dispersion than with the stellar mass measured 
within the fiber in the sense that the dispersion with the 
total mass is smaller that that of stellar mass in the fiber 
and Al/disp in Col. (4) is larger than that in Col. (7). 

In summary, only from the results of this paper, it is 
difficult to draw a robust conclusion whether the stellar mass 
or the total mass is the more fundamental parameter for 
the elliptical galaxies. But our results do suggest that the 
star formation history of nearby elliptical galaxy is primarily 
determined by the mass of the system. 

Our analysis of spectral absorption features has shown 
that both optical and optical-infrared colours are sensitive to 
variations in total stellar metallicity and, to a lesser extent, 
element abundance ratios. Age has a more significant effect 
on the optical colours than on the optical-infrared colours 
of ellipticals. The fact that optical-infrared colours are as 
strongly correlated with velocity dispersion as the optical 
colours argues for the fact that metallicity (and possibly el- 
ement abundance ratios) is the primary driver of the colour 
trends among elliptical galaxies. This is in agreement with 
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Figure 8. Left: g — r colour, Mg2/(Fe) and H7A index strengths 
against total mass Mtot = (2a meas ) 2 Rso,z /G . The black solid 
line in each panel represents the linear relation fitted for the sam- 
ple as a whole, while the coloured, dot- dashed lines show the re- 
lations for galaxies in four bins of stellar mass ( cyan: log M* < 
10.5, blue: 10.5 < logM, < 10.8, red: 10.8 < logAf* < 11.2, 
magenta: logAf* 11.2). Right: g — r colour, Mg2/(Fe) and 
H"YA index strengths against stellar mass. The coloured, dot- 
dashed lines show here the linear relations for galaxies in four 
bins of total mass (cyan: log M tot in j < 10.5, blue: 10.5 ^ 
log Mtot,inf < 10.8, red: 10.8 < log Mtot,inf < H-2, magenta: 
log M to t,inf ^ 11.2). In each panel the rms scatter about the 
relation is also indicated. 



the results bv lGallazzi et all i2005t) . who show that in con- 
trast to metallicity, the luminosity-weighted ages of early- 
type galaxies in SDSS are only weakly dependent on mass. 
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Table 4. Correlations between colours, spectral indices and the mass inferred from the measured velocity dispersion M totiri f = 
(2(T meas ) 2 -R50, z/G and stellar mass contained in the fiber M^jiber- Columns 2,5 give the slope of the relation; columns 3,6 give the 
dispersions about the relation; columns 4,7 give the change in colours (or index strength) along the fitted relation over the interval in 
structural parameter that contains 90 percent of the points, normalized to the dispersion about the fitted relation. 
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